arXiv: 1508.0553Ivl [math-ph] 22 Aug 2015 


Generalizing Bagarello’s operator approach to solve a class of partial 

differential equations 

Jean Ghislain Compaore, Villevo Adanhounme, and Mahouton Norbert Hounkonnou 


Abstract. The non-commutative strategy developed by Bagarello (see Int. Jour, of Theo¬ 
retical Physics., 43 , issue 12 (2004), p. 2371 - 2394) for the analysis of systems of ordinary 
differential equations (ODEs) is extended to a class of partial differential equations (PDEs), 
namely evolution equations and Navier-Stokes equations. Systems of PDEs are solved using an 
unbounded self-adjoint, densely defined, Hamiltonian operator and a recursion relation which 
provides a multiple commutator and a power series solution. Numerous examples are given in 
this work. 
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1. Introduction 

Evolution equations and linearized Navier-Stokes equations are extremely important in Math¬ 
ematics and Physics, especially for applications in engineering as they are often used to obtain 
vital information on flow phenomena such as the velocity distribution, the flow pattern and pres¬ 
sure losses. However, it is a well known fact that exact solutions of evolution equations are rather 
difficult to obtain. Equally, till date, very few exact solutions of the Navier-Stokes equations are 
known. This is largely due to the complexity of the system of differential equations involved. 

In the absence of a general solution, it is often convenient to experiment with models so 
as to obtain information on the flow phenomena. Chengri Jin and Mingzhu Liu used the 
Adomian decomposition method to solve a kind of evolution equation. Lu applied the “back¬ 
ground” method of the arbitrary Prandtl number problem to the Navier-Stokes equations in order 
to derive a scaling lower bound on the space-time averaged temperature of the layer along with 
an explicit prefactor and then used a multiple boundary layer asymptotic theory to sharpen the 
estimate, thus increasing the prefactor in the lower bound by another factor. The time-dependent 
discretized Navier-Stokes equations with Chorin’s projection scheme using the continuous Galerkin 
finite element method in time and space were studied by Kharrat and Mghazli[^. More recently, 
the existence of global regular solutions to Navier-Stokes equations was proven by Zadrzynska and 
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Zajaczkowski m while the exact solutions of rotating shallow water were obtained by Chesnokov 
[4] using group analysis. 

In this paper, we investigate the analytical solutions to two categories of PDE using Bagarello’s 
approach. The first category comprises evolution equations which arise in many modeling problems 
in physics and hydrodynamics. Well known examples include the KdV, MKdV, BBM, Burgers, 
regularized long wave (RLW) and KdV-KSV equations, which are all special cases of the general 
equation in [^. The second category is made up of Navier-Stokes equations and heat propagation 
equations. Equations in this category are used to describe phenomena in continuum mechanics 
and are derived from the mathematical model of the diffusion of particles in a medium. 

The rest of the paper is organized as follows. In Section 2, we present some details of the models 
we analyze. In Section 3, we obtain the analytical solutions to evolution equations, linearized 
Navier-Stokes equations and the pressure using a non-commutative approach. The validity of the 
technique is verihed through illustrative examples. This is followed by the conclusion. 


2. Analytical solution 

Our aim in this section is to solve all the types of equations enumerated above. 


2.1. Evolution equations. We consider the general evolution equation of the form 


du 

'dt 


M 


= E 


gn 


N 




E 


m—0 m—0 

subjected to the initial condition 


gm^k+l 

‘ g^m ^ dx'-dt 


+ f(x, t), — OO < X < -|-oo, t > 0 


i(x,t)\t=o = h{x) 


( 1 ) 

( 2 ) 


where / and h are given functions, a™, bm and c are real constants, k and i are positive integers, 
M and N are nonnegative integers. 

We suppose that the solution of the previous Cauchy problem exists and is unique. Then 
we solve this problem by Bagarello’s approach [^, based on a suggestion drawn from quantum 
mechanics. Indeed, given a quantum mechanical system S and the related set of observables Os, 
that is, the set of all the self-adjoint bounded (or more often unbounded) operators related to S, 
the evolution of any observable Y G Os satisfies the Heisenberg equation of motion (HOEM) 

^V(t,x) = i[ff,Y(t,x)]. (3) 

Here [A,i?] := AB — BA is the commutator between A,Bg Os, H is assumed to be a densely 
defined self-adjoint Hamiltonian operator of the system acting on some Hilbert space T-L. The 
initial condition Y^ = h{x) is also considered as an operator acting on the Hilbert space T-L. 

The following statement holds: 


Theorem 1. A formal solution of the HOEM is Y{t,x) = where Y^ is the initial 

condition ofY(t,x) and H does not depend explicitly on time. Furthermore if H is bounded, we 

get Y(t,x) = \H,Y^]k where [A,B]k is the multiple commutator defined recursively as: 

[A,B]o = B-,[A,B]k = [A, [A,B]k-i]. 


Proof. Indeed, we have 

dt' 
dt' 

m' 


-Y{t,x) 

-Y{t,x) 

-Y{t,x) 


= 

do 

= iHe^ 


'y°e- 


+ E‘ 








i[H,Y{t,x)\. 
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If H is bounded, Y{t,x) = can be expanded as 

k=0 

□ 


In order to apply Bagarello’s technique to the initial value problem we need the following 
Lemma 1. Assume 


n>0 


Then, we have 




{itY' 


n>0 


where < = XlLo ^nwlwl_\, p>2, = 


d 


Proof. Using the product of series 


2-^ 2-^ 2—j 2—d p\(n - 

j>0 fc>0 n>0 k=0 


k A n—k 


B^A 


we can write 

u^{x, t) 




j>0 k>0 


“ k\{n — k)\ 

[H,Y\ 


k\ 




n>0 k—0 


n>0 


Yx,t) = 


YY\jfixy'‘] 


j>0 k>0 


(*^) „„2 


fc! 


-Wk 




n>0 k—0 


n>0 


where = J2l=o and = J2l=o Cn^lK-k- 

By induction we get 


\x,t) = wl 


n>0 


n\ 


= p>2. 


fc=0 


In fact 


yielding 


,p+i 


{x,t) = u{x,t)uP{x,t) = 


1 


^ t!(„ 

n>0fc=0 ^ ’ 


/P+U 


{x,t) = 




n\ 

n>0 k—0 


□ 


Now we state the first of the main results through the following: 
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(it)^ 

Theorem 2. Ifu{x,t) = X]n>o —® solution of the problem dip, dii, for f{x,t) =0, 


then must satisfy the following recursive relation: 


Wn 




= Y° 

M 


'^n+1 


D I 


d” 


N 


J2 


5” 


m=0 


9a;" 


Vn > 0, 


ax* ^ "‘ax™ 

m=0 

w/iere = [H, y0]„, < = XlLo <^nwlwl~_\. 

Proof. Setting the expansions of u{x,t) and u^~^^{x,t) in the equation (P), we get 


( 4 ) 

( 5 ) 


E 

n>0 


{ity 


IW. 


n+1 


. 9* 1 


M 


-E 


on 


N 


m—0 


9x* 




9” 


which allows the following recursive relation 


9* 


= ro 

M 


m—0 


N 


dx' 


-w. 


fc+i 


= 0 , 


iw, 


n+1 


• ^ 1 

= E 


m—0 


'9x' 


7W„ 


E^” 

m—0 


dx' 


-w: 


k+l 


□ 


Remark. One can hnd wf, n = 1,2,3,... as particular solutions of the relation ([5|). 


2.2. Diffusion equations. Processes of heat propagation and of diffusion of particles in a 
medium are described by the general diffusion equation [1]: 


Ou 

p— - VfpVu) + qu = F{x, y, z, t), 


( 6 ) 


where p is the coefficient of porosity of the medium and p and q characterize its properties. In 
examining heat propagation, by u(x, y, z, t) we denote the temperature of the medium at the point 
(x, y, z) at time t. Assuming that the medium is isotropic, p(x, y, z), c(x, y, z) and fc(x, y, z) denote 
its density, specific heat capacity and heat conduction coefficient, respectively, while F(x,y,z,f) 
represents the intensity of heat sources. The process of heat propagation is described by a function 
which satisfies an equation of the type: 


( 7 ) 


( 8 ) 


Ou 

cp-^ - VfkVu) = T(x, y, z, f). 

If the medium is homogeneous, i.e. c, p and k are constant, equation 0 has the form: 

9 

—u = a'^Au + f{x,y,z,t), 

k F 

where = —, / = — and the last equation is the heat conductivity equation. For a complete 
cp cp 

description of the sources of heat propagation, it is necessary to specify the initial distribution 
of temperature u in the medium (initial conditions) and the conditions at the boundary of this 
medium (boundary conditions). 

Suppose that the equation ([5|) is subjected to the initial condition 

u\t=o = uo{x,y,z). (9) 

For given functions uq and /, one can solve the initial value problem using Bagarello’s approach. 
Here we formulate the second result: 


Theorem 3. Assume that uq G C°°(R^), the series 

+ 00 rfc 

E T7^^^o(a:,y,^:), <5 > 0 


/c=0 
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and those obtained by differentiating termwise up to second order converges uniformly on each 
bounded and closed subset C R^. Let 

+°° ( 2\k^k r 

ui{x, y, z,t)^'^ — - — A^uoix, y,z), 0<t< — 
k=0 

be a solution of the equation 

— = a^Au (10) 

satisfying the initial condition m where A is the Laplacian operator. 

If 

, . iitY 

U2{x, y, z,t) = ^ 
k>0 

is another solution of equation cni) satisfying the condition m, then we have 

Wk = {-ia‘^)^A'^uo{x,y,z), k>0, (11) 

where Wk = [H, Y'^]k and = Uo(x, y, z). 

Proof. Indeed, 


9m 1 


+°° (n'^Yik 

{x,y,z,t) = a'^Y^ — - — A'‘+^uo{x,y,z) 


/c=0 


9mi 

-^[x,y,z,t) 

9mi 

-^[x,y,z,t) 


= a^A 


(^2\k+k 


__ (^2\kfk 

Y — u — ^^'^oix,y,z) 


^ k^O 


= a'^Aui{x,y,z,t). 


Then, mi is a solution of such that Mi(x, y, z, t)|t=o = mo(x, y, z) . On the other hand, if M2 is 
a solution of m satisfying M 2 (x, y, z, t)|t=o = mo(x, y, z), we can write 


9m2 


+ 00 

{x,y,z,t) = Y 




j-k 


2^Au2 = ^ 


2 = 




which implies 


Wk+i = -ia^Awk- 


( 12 ) 


The last relation allows to write 

Mil = -ia^AY°, W2 = (-iaYl^^y°, W3 = {-ia^fA^Y°. 

By induction, we have Wk = (—By virtue of the uniqueness of the solution of Cauchy’s 
problem, we obtain ui{x, y, z, t) = U 2 (x, y, z, t) which implies 

i’^Wk = {a^)'"A'"uQ{x,y,z). 


Thus, Wk can be expressed as 

Wk = {-ia‘^)^A'^uo{x,y,z). 


That ends the proof. 


□ 
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3. Applications 

3.1. Example 1. We consider the RLW equation 
du 1 dv?' d^u 


dt 2 dx dx'^dt 
subjected to the initial condition 

u{x,t)\t=o = X. 

Using the recursive formula (O, we have 


, —oo < X < +0O, t > 0, 






Wq = X 

1 d 


2 dx 




(13) 

(14) 

(15) 

(16) 


k=0 


Finding the particular solutions of the previous recursive relation in the form = UnX, we obtain 
w\ = ix, W 2 = 2\i^x, wl = 4!i^x, ui| = 5!i^x. By induction we obtain = nli'^x. 


By virtue of the expansion of u{x, t), we get 

{it) 




n>0 


n>0 n>0 


-nil X 


+ t 


(17) 


We can easily see that (ITTl) is the exact solution of the initial value problem. 
3.2. Example 2. We consider 

d^n. 

— OO < X < +00, t > 0 


du du ^ d^u 

^ ^dx'^dt' 


u{x,t)\t=e = e 
Using the recursive formula ([S|), we obtain 


Wn = e ^ 


wi, n = 0,1,2,3,.... 


o - 1 1 9 

Thus, determining the particular solutions of the recursive formula we have 


(18) 

(19) 

( 20 ) 
( 21 ) 


w, = le 


Wn = i e 


w\=i^e-^\ 


4„-a: 


W 4 = I e 


By induction we obtain = z"e By means of the expansion of u{x, t), we can write 

/ J.\ (*^)" -n -X {i^t)'^ (~^)” -X -t-x 


d 


( 22 ) 


n>0 


n>0 n>0 _ 

as the exact solution of the initial value problem. 

3.3. Example 3. Find the solution of the initial value problem 

du ^d'^u d^u 

^ dx^dt 

u{x,t)\t=o = sin a;. 

By virtue of the recursive formula ([5]), we get 


.5^ 1 

OX^ 


Wq = Sin X 
= 2-®' 


n+1 ^ dx'^ 




(23) 

(24) 

(25) 

(26) 


Thus, w\ = isina;; w\ = i^sina;; w\ = i^sinx; w} = i^'^sinx; wl = i^sinx. By induction, we 
obtain wl, = sin x and 


*(a;,<) = Y 




I smx 


n>0 


n\ 


= E 


{iHY 


• smx 


n>0 


n\ 


= E 


{-^T . -t . 

-:— Sin X = e sin x 


n>0 


(27) 
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as the exact solution of the initial value problem. 


3.4. Example 4. Consider a homogeneous ball with constant radius R. The temperature T 
of this ball depending on the variable radius r and the time t is governed by the heat conduction 
equation [9] 


dT _ ^/d^T 2dT\ 


dt 


Var2 


dr 7 


subject to the initial condition 


and to boundary conditions 




T(r,t)|i=o = To(r) 


(r,t) + hT{rR)^\r=R = 0. 


In order to solve the latter problem, we undertake the change of variables V = rT and thus, 


dV 


,d^V 


= a 


dt dr^ 

V{r,t)\t=o = rTo{r) 

V{r,t)\r=o = 0; 


\r=R = 0. 


By virtue of the recursive relation (El), we have 


2 \k 


j2fc 


Wk = 

and the solution can be expressed in the form 

(it) 


dj-2k 


[rToir)] 


= E = E 




-[rToir)]. 


n>0 


n>0 


Going back to the initial variable, we get 


r^^ ^ {o?tY r / M 


r ^—' n\ 

n>0 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


3.5. Example 5. We consider the time-dependent Navier-Stokes problem in the primitive 
variables Ifillsl 


pAm -|- Vp 

= f 

in 

n X [0, T] 

(37) 

V.u 

= 0 

in 

n X [o,T] 

(38) 

u 

= 0 

on 

rx]o,T[ 

(39) 

U\t=0 

= Uo 

in 


(40) 


where fl is the bounded connected domain of {d = 2,3) with the Lipschitz continuous boundary 
T; the unknowns are the velocity u = {ui{x,y, z,t),U 2 {x,y, z,t),U 3 {x,y, z,t)) and the pressure 
p = p{x,y,z,t)-, the data stand for / = {fi{x,y,z,t),f 2 (x,y,z,t),f^{x^y,z,t)) which represents 
the force on the body and uq = ito(a;, y, z), the initial condition; v is the kinematic viscosity 
assumed to be a positive constant. 

Using Bagarello’s approach, we provide analytical solutions to the linearized Navier-Stokes 
equations (1571) satisfying the initial and boundary conditions (I551) - (H(I1) and the pressure. 

We can state the main result as: 










8 JEAN GHISLAIN COMPAORE, VILLEVO ADANHOUNME, AND MAHOUTON NORBERT HOUNKONNOU 


Theorem 4. The velocity u and the pressure p of the incompressible fluid flow described by 
the Navier-Stokes equations are defined by 


where 


u{x,y,zfi) = -V X A '^fi{x,y,zfi)+ Vtp[x,y,z,t) in rj x [0, T] 
p{x,y,z,t) = po + 2 / 0 , 20 ,io) - V.(A-V)(a;o, 2 / 0 , ^o,io) 

-^(a;, 2/, 2,0 + V.(A~ V)(a;, 2/, 2, <), 


ip{x,y,zfi) = '^—^lsT[Vxuo{x,y,z)\+(j){x,y,zfi) 


n>0 
pT 




A ^v{x,y,zfi) ~ J J 

w = ix,y,z), 

(j) a particular solution of 

^-lyAtfj = Vx/ in f2x[0,T], 


v{^, t)dfdr 


(41) 

(42) 

(43) 

(44) 


(45) 


(fi is a given harmonic potential and its derivatives are decreasing functions when \/-V y^ z^ 
tends to +oo. 

Proof. Indeed, applying the operator V to equation dSZl) and setting if = V x u, we have 

dif 


dt 


-uAif = Vx/ in rJx[0,T] 
fi\t=o = V X itQ in 


(46) 

(47) 


Setting V X / = 0, substituting if = 


{ity 

^0 n! 


the previous theorem, we have 

[7^,r°]„ = (-^^.)"A"T^ ^(a;,y,z,t) = ^i^A"r°, = V x uq- (48) 


[H,YX in the equation (1461) and by virtue of 

i^ty 


n>0 


Therefore, we can express the solution 'll; to the problem (j46|l and (|Tf|l by 
yix, 2/, i) = T - j- A”F° + (f{x, y, z, t), 

Z—^ 77,1 


(49) 


n>0 


where </> is a particular solution of equation (H51) . 

Applying the operator V to equation if = V x u and taking account of eauation (l38l) . we obtain 
the equation with unknown u 


Au = —V X if 

which allows to get the solution u in the form [8] 

u{x, y, z, t) = -V X A~^if{x, y, z, f) + V(p{x, y, z, t), 
where the inverse operator A“^ of A is given by 


A 


{x,y,z,t) = ^(dTrr) exp | '^v{^,t)dfdT, uj = {x,y,z), 


(50) 

(51) 

(52) 
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ip a given harmonic potential of the fluid and its derivatives are decreasing when + y'^ + z'^ 
approaches +oo. Then we derive the pressure p from eauation (l37D as 
^ du 

ot 
d 

Vp = / — ^( — Vx + Vp) — vV X ij; 


/-|(V„) + A->(vx|:)-.Vx^ 

d 


V X 


dif) 


— I 


/-^(Vp) + A-i Vx Vx/ 
at L 


= V 


dp 

'm 


+ V.{A-\f) 


which reduces to 


p=—^+V.(A ^/l + C, C = const. 
at 


(53) 


For measured stagnation pressure p{x,y,z,t) = po at the point {xo,yo, Zo^to), the pressure p 
exerted on the fluid flow is defined by 

p{x,y,z,t) = po + -^{xo,yo,ZQ,tQ) -\/.{A~^f){xo,yo,zo,to) 

-^{x,y,z,t) + V.{A~^f){x,y,z,t). 

That ends the proof. 

Without loss of generality, we set 

Vxuo{x,y,z) = {cosy cos z,sm{x — y — z),e^'^'^^^), 

Vxf{x,y,z,t) = {tcosx,e*,tzsmx) 

and thus obtain 


□ 


cosy cosz + ly ^{—1 + iyt + e ‘"*)cos: 
tjj{x, y, z,t) = I sin(a; — y — z) + e* — I 

^3,^t+x+y+z _P _P e-‘'*)z sinx 


^3ut+x+y+z ^-3ut -y-z) 

V X tjj{x,y, z,t) = ( —cosj/sinz — 1 + z/t + cos: 

g-2!/t y (,Qg ^ _ y _ 2 ;'j 


A ^{V X ip{x,y,z,t)) = j ^(47^'^) ^'^^expj '^{'^ x'il;{^,t))d^dT. 


Choosing the harmonic potential ip in the form 

t 


p{x,y,z,t) = 


\/x2 +y‘^ + z^ 
we can express the velocity u and the pressure p as 


{x,y,z) ^ (0,0,0), 


(54) 

(55) 

(56) 

(57) 


u{x,y,z,t) = - 


J J (47rr) ^/^exp |—x tp{^,t))d^d7 


P{x,y,z,t) = pq + 


-t[x{x'^+ y'^ + z'^) ,y{x^ + y^ + z^) z(a;^ + + z^) 

1 1 


\/xl+yl + zl ^/sF+Tp+Z^ 
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4. Concluding remarks 

The non-commutative strategy developed by Bagarello for the analysis of systems of ordinary 
differential equations (ODEs) has been extended to a class of partial differential equations. Sys¬ 
tems of PDEs have been solved using an unbounded self-adjoint, densely defined, Hamiltonian 
operator and a recursion relation which provides a multiple commutator and a power series solu¬ 
tion. Numerous examples like the equations of diffusion, RLW, heat conduction and Navier-Stokes 
equations have been given in this work proving the efficiency of considered method also for the 
analysis of relevant PDEs widely used in mathematical physics. In that sense, this work can be 
considered as a good supplement to Bagarello’s approach. 
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